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List  of  Variables 


Numerous  variables  are  used  throughout  this  report.  These  variables  are 
listed  here  for  convenience.  Figures  1  and  2  show  the  measurements  that  define 
the  variables  involving  the  leaf  mass  and  plant  dimensions  for  submerged  and 
emergent  (unsubmerged  or  partially  submerged)  flow  conditions. 

The  variables  are  defined  as  follows: 

A  Cross  sectional  flow  area,  ft2  or  m2 

Ai  Frontal  area  of  an  individual  plant  blocking  flow,  approximated  by 
the  equivalent  rectangular  area  of  blockage  H’  by  We,  ft2  or  m2 

Ai *  Net  submerged  frontal  area  of  a  partially  submerged  plant,  ft2  or  m2 

As  Total  cross-sectional  area  of  all  of  the  stem(s)  of  an  individual  plant, 
measured  at  H/4,  ft2  or  m2 

b  width  of  channel  flume,  ft  or  m 

C  Chezy  resistance  coefficient,  fd/i/s  or  m  Ws 

CD  Drag  coefficient  of  vegetation,  dimensionless 

Ds  Stem  diameter,  measured  at  a  height  of  H/4,  ft  or  m 

Es  Modulus  of  plant  stiffness,  lbf/ft2  or  N/m2 

/  Darcy-Weisbach  friction  factor,  dimensionless 

fb  friction  factor  for  the  bed  and  plants,  dimensionless 

fw  friction  factor  for  the  walls,  dimensionless 

F45  The  horizontal  force  necessary  to  bend  a  plant  stem  45  deg,  lbf  or  N 

FD  Drag  force,  lbf  or  N 

Fr  Froude  number,  dimensionless 


g  Acceleration  due  to  gravity  =  3217  ft/s2  or  9806  m/s2 

H  Average  undeflected  plant  height,  ft  or  m 

If  Undeflected  height  of  the  leaf  mass  of  a  plant,  ft  or  m 

H*  Undeflected  height  of  leaf  mass  that  is  below  water  surface  for  a 
partially  submerged  plant,  ft  or  m  (See  Figure  2) 

I  Second  moment  of  inertia  of  cross  section  of  plant  stem,  ft4  or  m4 

Kn  Units  conversion  factor  for  Manning’s  equation,  14861  ftl/3/sec  or 
10  ml/3/sec 

L  Channel  reach  length,  ft  or  m 
M  Relative  plant  density,  number  of  plants  per  ft2  or  m2 
n  Total  Manning’s  roughness  coefficient,  including  sidewall  roughness 
rib  Manning's  resistance  coefficient  for  vegetation  and  channel  bed 
ttyeg  Manning's  resistance  coefficient  for  vegetation 
na  Manning’ s  resistance  coefficient  for  the  bed 
P  Wetted  perimeter,  ft  or  m 

Re  Reynolds  number,  R*  =  V  Rh  A 

Rh  Hydraulic  radius,  Rh  =  flow  area  /  wetted  perimeter,  ft  or  m 
Rb  Hydraulic  radius  for  the  bed  and  plants,  ft  or  m 
Rw  Hydraulic  radius  for  the  walls,  ft  or  m 

S  Bed  or  energy  slope,  dimensionless 

S0  Bed  slope ,  dimensionless 

Sf  Energy  slope,  dimensionless 

V  Mean  channel  velocity,  ft/s  or  m/s 

Vp  Local  plant  approach  velocity  in  front  of  the  leaf  mass,  ft/s  or  m/s 
V*  Shear  velocity,  V*  =  (g  Rh  S)V2 ,  ft/s  or  m/s 
V*/V  Resistance  coefficient,  dimensionless 


Y0  Flow  depth,  ft  or  m 

We  Equivalent  average  plant  width,  We  =  A,/H\  ftorm 
dy/dx  Unit  change  in  slope  of  the  water  surface 

y  Specific  weight  of  water,  lbf/ft3  or  N/m3 

v  Fluid  dynamic  viscosity,  ft2/s  or  m2/s 

p  Fluid  density,  slugs/ft3  (lbf-sec/ft)  or  kg/m3 

to  Shear  stress  on  channel  bottom,  (to  =  y  Rh  S),  lbf/ft2  or  N/m2 


1  Introduction 


Background 

An  important  consideration  for  determining  the  stage-discharge  relationship 
in  rivers  and  streams  is  the  effect  or  influence  of  vegetation  on  the  overall  head 
loss  along  a  channel  and  in  the  overbank.  Plants  in  the  floodplain  and  along  the 
banks  can  increase  or  even  decrease  the  effective  flow  resistance.  The  vegetation 
may  be  natural  or  it  may  have  been  planted  to  improve  aesthetics  or  habitat,  to 
prevent  erosion,  or  for  other  reasons. 

The  impetus  for  this  study  came  as  a  result  of  numerous  inquiries  from 
U.S.  Army  Corps  of  Engineer  District  offices  regarding  the  proper  hydraulic 
roughness  values  to  use  for  shrubs  and  other  aesthetically  and  environmentally 
desirable  plants.  The  District  offices  were  involved  in  the  evaluation  of  vegeta¬ 
tive  impacts  on  proposed  and  existing  channels  to  determine  flow  capacity  and 
water  surface  elevations.  Given  the  near  complete  lack  of  hydraulic  roughness 
values  for  shrubs  and  similar  vegetation,  the  accurate  estimation  of  channel 
capacity  and  water  surface  elevations  was  difficult  at  best.  The  work  described 
herein  was  carried  out  under  the  Flood  Control  Channels  and  Flood  Damage 
Reduction  Research  Programs  starting  in  1993  and  completed  in  1997.  It  was  a 
direct  result  of  District  requests  for  research  through  the  Flood  Control  Channels 
Field  Review  Group 

Previous  research  has  been  conducted  on  vegetation  such  as  grasses,  agricul¬ 
tural  crops,  and  on  the  rigid  blockage  of  cylindrical  tree  trunks.  However,  little 
had  been  studied  on  the  resistance  effects  of  plants  and  shrubs  that  are  either 
submerged  or  partially  submerged  by  turbulent  flows.  The  flexible  stems  and 
varying  shapes  of  plant  leaf  mass  greatly  complicate  the  understanding  of  this 
resistance.  The  deformation  of  plant  shape  with  flow  precludes  the  use  of  a 
constant  blockage  area  or  the  density  of  plant  frontal  area  in  predicting 
resistance. 

The  purpose  of  this  study  was  to  investigate  the  effect  of  vegetation, 
particularly  ground  cover  plants,  small  trees,  and  shrubs,  on  flow  resistance. 
Hydraulic  losses  and  drag  due  to  actual  plants  were  measured  at  the  Utah  Water 
Research  Laboratory  utilizing  a  large  wide  flume  and  a  smaller  sectional  flume. 
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Research  in  the  flume  resulted  in  the  collection  of  data  from  more  than 
220  experiments  with  20  different  plant  species.  Experiments  were  conducted 
with  both  homogeneous  and  mixed  plant  groupings.  Single-stem  and  multiple- 
stem  plants  were  included  in  the  plant  types  evaluated.  Plants  with  and  without 
leaves  were  evaluated.  Plant  density,  spacing,  and  size  were  varied  in  the 
experiments.  Plants  were  evaluated  over  a  range  of  velocities  and  depths. 

A  methodology  was  developed  from  the  laboratory  data  to  predict  head  loss 
and  resistance  coefficients  as  a  function  of  slope  and  depth.  Input  data  for  the 
methodology  can  be  collected  from  the  field  or  estimated  plant  characteristics 
may  be  used. 


Resistance  Coefficients 

Resistance  to  flow  is  typically  characterized  by  a  roughness  coefficient.  The 
most  commonly  used  equation  for  flow  resistance  is  the  Manning's  equation: 

K  V  V 

V  =  —R(3  S/2  (1) 

n 

where 

V  =  mean  velocity  of  flow 

Rh  =  hydraulic  radius 

5  =  slope  of  the  energy  grade 

n  =  Manning’s  resistance  coefficient 

Ka  =  unit  correction  factor  of  1 .0  for  SI  units  and  1 .486  for  non-SI  units 

Although  Manning’s  equation  is  used  extensively  for  calculating  flow 
resistance,  Manning  himself  did  not  recommend  it  for  use,  because  his  research 
found  that  n  was  not  constant  but  varied  with  velocity  and  depth. 

The  ratio  of  shear  velocity  to  mean  velocity,  V*/V,  is  another  form  of 
resistance  coefficient.  Keulegan  (1938)  used  it  to  calculate  average  velocity 
based  on  the  theoretical  vertical  velocity  profile.  The  ratio  of  shear  velocity  to 
average  velocity  may  be  thought  of  as  the  ratio  of  shear  stress  to  inertial  force  as 
indicated  in  Equation  2. 


V*.  _  jgR^S 

V  V 


(2) 


V*  is  the  shear  velocity,  g  is  the  acceleration  due  to  gravity,  T0is  the  shear  stress, 
and  p  is  the  density  of  water. 
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There  are  other  resistance  coefficients  in  use  including  the  Darcy-Weisbach 
friction  factor,  f,  and  the  Chezy  C.  These  can  be  converted  easily  to  Manning’s 
n  as  shown  in  Equation  3. 


(3) 


Note  that  the  Chezy  coefficient  is  not  dimensionless,  and  will  vary  with 
units.  The  Manning  coefficient  is  dimensionless  only  in  the  sense  that  Kn  will 
make  it  dimensionless  when  units  of  feet  or  meters  are  used  in  Equation  3. 


In  this  study,  resistance  equations  were  developed  for  the  shear  velocity  to 
average  velocity  ratio  because  it  is  dimensionless  and  has  a  sound  theoretical 
basis,  and  for  the  Manning’s  coefficient  because  its  use  is  widespread. 

The  Manning’s  resistance  coefficient  for  vegetation  is  calculated  in 
conformity  with  the  Cowan  (1956)  method  for  additive  resistance.  This  method 
consists  of  additions  to  roughness  for  various  surface  irregularities  and 
vegetation.  The  equation  that  describes  the  method  is: 

n  =  (n0  +nj  +n2  +n3  +n4  )m  (4) 


where  n„  is  the  base  value  for  a  straight,  uniform,  smooth  channel  in  natural 
materials;  ni  is  an  additive  value  to  account  for  surface  irregularities;  n2  is  added 
to  account  for  variations  in  the  channel  geometry  along  the  reach;  n5  is  an 
additive  value  to  account  for  obstructions;  a*  accounts  for  vegetation;  and  m  is  a 
correction  factor  for  meandering  or  sinuosity  of  the  channel.  The  114  coefficient 
used  in  Cowan’s  method  is  based  on  the  net  effect  of  vegetation. 


Many  published  values  of  Manning’s  roughness  coefficients  related  to 
vegetated  surfaces  include  the  base  resistance,  no,  as  a  part  of  the  reported 
vegetation  resistance.  This  is  the  convention  followed  in  this  report.  Thus, 
roughness  coefficients  reported  herein  include  the  effects  of  both  the  bed  and  the 
vegetation.  In  Cowan  notion  this  would  be  expressed  as  n=n0+n4. 

The  resistance  values  can  be  composited  into  channel  averages  using  several 
methods.  One  of  the  methods  is  Lotter’s  1933  method  shown  in  Equation  5.  The 
Lotter  method,  presented  by  Chow  (1959),  uses  vertical  bisecting  lines  to 
subdivide  the  channel  cross  section  into  subareas  for  the  calculation  of  flow. 

This  method  assumes  that  the  total  flow  is  the  sum  of  the  flows  in  the  separate 
subareas.  The  equivalent  resistance  thus  developed  accounts  for  the  variability  in 
resistance  across  the  channel. 


n  = 


PRhA 


i= 1 


Pi  Rh-A  A 


n,. 

v  J 

P  is  the  wetted  perimeter.  N  is  the  number  of  i  subsections. 


(5) 
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Substitution  of  Equation  5  into  Equation  3  yields  a  dimensionless  composit¬ 
ing  equation  for  V*/V  given  in  Equation  6. 


V* 

V 


Pr/2 


(6) 


Resistance  Equation  for  Large  Woody  Vegetation 

Usually,  vegetation  on  the  flood  plains  is  larger  than  that  found  in  the  main 
channel.  This  vegetation  has  a  major  influence  on  flow  depth  and  resistance 
during  overbank  flooding.  Petryk  and  Bosmajian  (1975)  proposed  a  method  to 
calculate  flow  resistance  based  on  the  drag  forces  created  by  the  larger  plants  and 
trees  that  constitute  much  of  the  resistance  on  the  flood  plains.  They  derived 
Equation  7  for  Manning's  n  by  summing  the  forces  in  the  longitudinal  direction. 
The  forces  include  pressure  forces,  the  gravitational  force,  shear  forces,  and  the 
drag  forces. 


n 


1.0  + 


Cv  ^  A 

2gAL 


vO’ 

n„ 


(7) 


Here  CD  is  the  effective  drag  coefficient  for  the  vegetation  in  the  direction  of  the 
flow,  A  is  the  cross-sectional  area  of  the  flow,  in  ft2,  XAj  is  the  total  frontal  area 
of  vegetation  blocking  the  flow  in  the  reach,  in  ft2,  L  is  the  length  of  the  channel 
reach  being  considered,  in  ft.  The  expression  CD  EAj/(AL)  represents  the 
vegetation  blockage,  or  the  density  of  vegetation  in  the  floodplain.  This 
expression  must  be  either  directly  or  indirectly  measured.  The  total  boundary 
roughness,  no,  excludes  the  additive  resistance,  114  for  other  types  of  vegetation 
such  as  shrubs. 


There  are  several  limitations  to  using  Petryk  and  Bosmajian's  equation.  The 
channel  velocity  must  be  small  enough  to  prevent  bending  or  distortion  of  the 
vegetation,  and  large  variations  in  velocity  cannot  occur  across  the  channel. 
Vegetation  such  as  grasses  and  shrubs  are  then  excluded.  Vegetation  must  also 
be  distributed  relatively  uniformly  in  the  lateral  direction.  Finally,  according  to 
Petryk  and  Bosmajian,  the  flow  depth  must  be  less  than  or  equal  to  the  maximum 
vegetation  height.  During  flooding,  the  velocities  over  the  floodplains  can  be 
relatively  high  and  large  degrees  of  bending  and  distortion  of  vegetation  often 
occur.  Vegetation  types  and  densities  can  also  vary  widely  across  a  floodplain, 
and  water  depths  often  submerge  vegetation.  However,  when  tree  trunks 
dominate  sections  of  a  floodplain,  this  method  can  be  used  to  predict  the  total 
resistance  coefficient,  n. 
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Laboratory  Setup  and 
Procedures 


Experimental  Plants 

Two  flumes  at  the  Utah  Water  Research  Laboratory  were  used  for  laboratory 
experiments  during  this  study.  The  large  flume  (Figure  3),  2.44  m  (8  ft)  wide  by 
1.82  m  (6  ft)  deep  by  152.4  m  (500  ft)  long,  was  used  to  measure  in  situ  flow 
resistance  and  drag  force  for  groups  of  uniform  sized  plants  and  groups  of  mixed 
plants  with  varying  plant  density,  sizes,  and  shapes.  A  sectional  flume,  0.91  m 
(3  ft)  wide  by  0.91  m  (3  ft)  deep,  was  used  to  measure  drag  force  of  individual 
plants.  Thirteen  different  plant  types  were  evaluated  in  the  large  laboratory 
flume  and  10  plant  types  in  the  sectional  flume.  Six  combinations  of  plants 
typical  of  different  ecosystems  were  also  studied  in  the  large  flume.  In  total, 

21  different  plant  types  were  evaluated  in  the  two  flumes.  The  plants  and  their 
characteristics  are  listed  in  Table  1.  Field  measurements  of  plant  stiffness 
(Freeman  1997)  for  four  plant  types  are  also  listed  in  Table  1.  All  plants 
evaluated  were  broadleaf  deciduous  vegetation  commonly  found  in  floodplain 
and  riparian  zones. 


DOWNSTREAM  TEST  ARTIFICIAL  MFFCSERS  INLET 

GATE  PLANTS  ROUGHNESS  AND 

BLOCKS  CONTROL 

VALVE 


Figure  3.  Layout  of  large  flume  for  plant  roughness  experiments  (To  convert  feet  to 
meters,  multiply  by  0.3048) 
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The  plants  evaluated  in  the  large  flume  were  placed  in  staggered  rows  along 
the  22.9-m  (75-ft)  length  of  the  experimentation  section  (Figure  3).  The  spacing 
selected  for  the  plants  was  based  on  typical  plant  spacing  found  in  floodplains. 
The  plant  density,  M,  was  calculated  as  the  number  of  plants  per  unit  area.  The 
plants  evaluated  in  the  small  flume  were  placed  in  a  single  row  of  four  to 
five  plants  along  the  center  line  of  the  flume.  A  single  plant  was  instrumented 
for  measuring  drag  force  in  both  flumes.  The  instrumented  plant  in  the  larger 
flume  was  located  in  the  center  of  the  22.9-m  (75  ft)  by  2.44-m  (8-ft)  experimen¬ 
tation  section.  The  plant  selected  for  measurement  in  the  small  flume  was  the 
downstream  plant,  with  four  plants  located  upstream.  The  experimental  setup  for 
the  small  flume  allowed  for  a  more  accurate  measurement  of  plant  approach 
velocity  (VP)  and  drag  force  (FD).  Roots  had  to  be  removed  from  all  the  plants 
used  in  the  small  flume,  but  only  the  plant  used  to  measure  drag  force  in  the  large 
flume  required  root  removal.  All  other  plants  in  the  large  flume  were  placed 
intact,  with  root  structure  and  original  soil,  into  a  20.3-cm  (8-in)  deep 
experimental  bed. 

The  range  of  variables  measured  in  the  large  flume  were: 

a.  Flow  depths  from  0.4  to  1.4  m  (1.3  to  4.7  ft). 

b.  Average  flow  velocities  from  0.15  to  1.1  m/s  (0.5  to  3.6  ft/s). 

c.  Measured  resistance  V»/V  from  0.13  to  0.45  and  n  from  0.04  to  0.14. 

d.  Plant  heights  from  0.20  to  1.52  m  (0.66  to  5  ft). 

e.  Plant  widths  from  0.076  to  0.91  m  (0.25  to  3  ft). 

f  Plant  densities  from  0.53  to  1 3  plants  /  m2  (0.05  to  1 .2  plants  /  ft2). 

g.  Plant  modulus  of  stiffness  from  5.3  x  107  to  4.8  x  109  N/m2  (1.1  x  106to 

1.0  x  108lbf/ft2). 

h.  Reynolds  numbers  from  1.4  x  105  to  1.6  x  106. 


Large  Flume  Experimental  Setup 

The  large  flume  experimentation  section  had  a  bed  that  would  accept  plants 
with  their  root  systems  intact.  In  the  first  phase  of  the  study,  the  bed  consisted  of 
a  layer  of  gravel,  to  assist  in  drainage,  covered  with  a  cap  of  compacted  clay. 

This  material  supported  the  plants  and  was  intended  to  prevent  them  from 
washing  downstream.  In  the  second  phase  of  the  study,  the  compacted  clay  was 
replaced  with  a  gravel  bed  and  a  mortar  cap.  The  mortar  cap  greatly  facilitated 
the  changing  of  plants  and  experimental  setups.  The  mortar  cap  had  158  3.8-L 
(1-gal)  plant  containers  in  the  bed  in  staggered  rows  of  four  and  five  plants  per 
row.  When  particular  containers  were  not  in  use  they  were  capped  flush  with  the 
top  of  the  mortar  cap  to  prevent  the  introduction  of  additional  roughness. 
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Upstream  and  downstream  of  the  experimentation  section  the  flume 
contained  a  section  of  roughened  bed.  The  roughness  elements  consisted  of 
cinder  blocks  that  were  adjusted  until  they  produced  a  fully-developed  turbulent 
velocity  distribution  upstream  and  downstream  of  the  experimentation  section. 

At  the  beginning  of  each  experiment,  at  the  downstream  end  of  the  clay  or 
mortar  bed,  stop  logs  were  inserted  to  allow  for  slow  filling  of  the  flume.  This 
was  done  to  protect  plants  during  filling.  As  discharge  was  slowly  increased  to 
the  desired  level,  stop  logs  were  removed.  Some  stop  logs  remained  during  the 
experiment  to  maintain  a  constant  velocity  profile  throughout  the 
experimentation  section.  At  the  downstream  end  of  the  flume,  91.4  m  (300  ft) 
downstream  of  the  experimental  section,  a  hydraulic  gate  was  used  to  control 
flow  depth. 

Water  from  the  river  adjacent  to  the  laboratory  entered  the  flume  from  a 
122-cm  (48-in)  pipe  50.3  m  (165  ft)  upstream  from  the  experimentation  section. 
Water  temperature  was  measured  and  found  to  be  10°C  (50°F)  for  all  experi¬ 
ments.  A  remote  controlled  butterfly  valve  in  the  122-cm  (48-in)  pipeline  was 
used  to  control  the  flow  of  water  into  the  flume.  A  Mapco  sonic  meter  was  used 
to  measure  the  flow  rate  in  the  inlet  pipe.  Downstream  from  the  inlet  pipe,  jet 
flow  was  dissipated  using  a  series  of  vertical  and  horizontal  distribution  vanes. 

A  wheeled  platform  that  moved  on  tracks  adjacent  to  the  flume  sides  was 
used  to  take  depth  and  velocity  measurements.  This  platform  was  positioned  at 
1.52-m  (5-ft)  intervals  along  the  length  of  the  experimental  section  to  facilitate 
measurements.  Water-surface  elevations  were  measured  with  the  help  of  a 
stationary  transit  and  a  measuring  rod  along  the  center  line  of  the  flume.  Flow 
velocities  were  taken  with  a  Marsh  McBimey  Model  201  portable  water  current 
meter. 

A  single  plant  centered  horizontally  in  the  flume  was  selected  to  measure 
drag  force.  An  average-sized  plant  was  selected  and  inserted  into  a  platform. 

The  platform  was  a  shallow  metal  box  with  ball  bearings  in  the  bottom  and  a 
metal  plate  resting  upon  the  ball  bearings  as  shown  in  Figure  4.  The  instru¬ 
mented  plant,  with  its  roots  removed,  was  attached  to  the  plate.  A  Vishay 
Instrument  Model  P-350  strain  indicator  was  then  attached  to  the  downstream 
end  of  the  plate  to  measure  the  drag  force  applied  to  the  plant  by  the  moving 
water  column.  This  drag  force  was  measured  as  a  compression  force.  During  the 
experiment  the  platform  was  covered  with  a  section  of  drain  cloth  to  prevent  soil 
from  interfering  with  the  ball  bearings  and  movement  of  the  plate.  The  platform 
was  also  covered  with  a  plastic  lid  to  reduce  friction  drag  on  the  platform.  The 
strain  gage  was  zeroed  at  the  start  of  each  series  of  experiments.  The  range  of 
the  strain  gage  was  0  to  44.5  N  (10  lbf).  The  sensitivity  of  the  strain  gage  was  45 
micro-cm  per  cm  per  N  (200  micro-inches  per  inch  per  pound).  Measurements 
were  taken  to  the  nearest  2.5  micro-cm  (1.0  micro-inch.) 
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Large  Flume  Operating  Procedures 

Measurements  were  made  of  plant  dimensions  and  plant  characteristics 
before  each  series  of  experiments.  Plant  height  and  width,  leaf  size,  and  stem 
height  were  measured.  The  number  of  branches  and  stems  were  counted.  The 
diameter  of  stems  and  branches  was  recorded,  and  bending  characteristics  were 
measured.  The  forces  required  to  bend  the  plant  45  deg  and  horizontal  at 
different  heights  along  the  stem  were  determined.  A  strain  gage  was  first 
attached  to  the  top  of  the  plant.  After  the  bending  forces  and  deflection  were 
determined  there,  the  gage  was  hooked  to  the  center  of  the  plant  and  the  bending 
forces  were  again  measured. 

During  the  first  phase  of  the  study,  prior  to  beginning  each  series  of  experi¬ 
ments,  the  bed  was  leveled  and  a  layer  of  topsoil  placed  and  compacted  on  top  of 
the  clay  bed.  The  mortar  cap  used  in  second  phase  of  the  study  did  not  require 
maintenance  and  leveling  for  each  series  of  runs.  The  plants  were  placed  in  the 
flume  just  prior  to  the  experiment  and  the  flume  was  slowly  filled  with  water, 
with  the  stop  logs  in  place  and  the  downstream  gate  closed.  With  the  flume 
filled  and  no  flow,  the  strain  gage  was  zeroed.  Flow  and  depth  were  controlled 
with  the  downstream  gate  and  the  122-cm  (48-in)  inlet  butterfly  valve.  Time  was 
allowed  for  the  flume  to  reach  equilibrium  before  measurements  were  taken  for 
each  run. 

Typically,  nine  runs  were  made  for  each  series  of  experiments.  The  first 
three  runs  were  made  at  high  depths,  with  the  flume  nearly  full,  and  at  three 
different  velocities.  The  next  three  runs  were  made  at  a  medium  depth,  and  the 
last  three  runs  were  made  at  a  low  depth.  The  plants  were  usually  submerged, 
even  at  low  depths,  because  the  flow  forces  were  adequate  to  bend  the  plants 
with  the  flow.  Some  runs  were  conducted  near  the  end  of  the  study  with  the 
larger  plants  partially  submerged,  to  aid  in  determining  a  relationship  for  partially 
submerged  plants. 
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Water-surface  elevation  was  measured  at  1.52-m  (5-ft)  intervals  along  the 
length  of  the  flume’s  experimentation  section.  At  the  midpoint  of  the  experimen¬ 
tal  section,  velocity  measurements  were  taken  at  different  depths  to  establish  a 
vertical  velocity  profile.  Velocity  measurements  were  also  taken  at  the  center  of 
the  leaf  mass  just  upstream  of  the  plant  used  to  measure  drag  force.  The  plant 
approach  velocity  was  measured  5.1  cm  (2  in)  upstream  of  the  instrumented  plant 
to  avoid  making  a  measurement  in  what  could  have  been  a  stagnation  region  at 
the  upstream  face  of  the  plant.  It  was  determined  early  in  the  study  that  velocity 
measurements  taken  in  the  plant  mass  and  at  the  upstream  face  of  the  plant  were 
inconsistent  because  of  the  interference  of  individual  leaves,  but  the  velocity 
measurements  did  show  that  there  was  still  substantial  velocity  and  flow  through 
the  plant  mass.  Drag  force  was  determined  from  the  strain  gage  measurements 
for  many,  but  not  all,  experimental  runs.  As  the  depths  and  velocities  were 
varied,  the  plants  and  bed  (for  Phase  I)  were  observed  through  the  view  window 
to  document  soil  movement,  plant  distortion,  and  plant  failure. 


Small  Sectional  Flume  Setup 

A  smaller  sectional  flume  was  used  to  study  the  drag  forces  developed  on 
single  plants.  Water  was  supplied  by  a  0.914-m-  (3-ft-)  wide  by  0.914-m-  (3-ft-) 
high  channel  running  perpendicular  to  the  flume  entrance.  A  baffle  was  placed  at 
the  entrance  of  the  flume  to  straighten  the  incoming  flow  and  a  Plexiglas 
observation  window  was  installed  in  the  side  of  the  flume. 

Since  the  bottom  of  the  flume  consisted  of  smooth  steel,  it  was  necessary  to 
devise  a  method  to  secure  the  plants  in  the  flume.  A  3.81-cm-  (lVi-in.-)  thick 
false  deck  was  constructed  of  smooth,  painted  plywood.  The  deck  was  bolted 
through  the  bottom  of  the  flume  and  sealed  with  silicon  caulk.  Several  2.54-cm 
(1.0-in)  holes  were  drilled  through  the  plywood  to  the  steel  bottom.  These  holes 
were  designed  to  hold  plants  in  a  layout  that  would  create  a  flow  regime  around 
the  plant  similar  to  the  flow  regime  of  the  plant  in  the  large  flume. 

To  attach  the  plants  to  the  flume  bottom,  a  beveled  rubber  grommet  and 
wide-flanged  washers  were  used.  The  roots  of  the  plants  were  cut  off  at  the  base 
of  the  stem,  and  the  plant  stem  inserted  through  the  washer  into  the  grommet. 

The  rubber  grommet  was  used  to  protect  the  base  of  the  stem  and  prevent 
breakage  of  the  stem.  Without  the  grommet,  the  plant  tended  to  break  where  the 
stem  contacted  the  surface  of  the  plywood  floor  when  the  plants  were  subjected 
to  high  velocities.  The  rubber  in  the  grommet  would  give  slightly,  however, 
allowing  the  plant  to  bend  a  small  amount  at  the  base  rather  than  shear  off  against 
the  sharp  edges  of  the  plywood  floor.  This  is  similar  to  the  conditions  that  the 
plant  experiences  in  the  field  with  soil  around  its  base.  The  wide  flanged 
washers  had  two  holes  that  allowed  the  grommet  to  be  attached  to  the  plywood 
floor  with  screws.  The  beveled  grommet  was  slightly  larger  than  the  holes  and 
when  the  screws  were  tightened,  the  washer  compressed  the  grommet  into  the 
hole,  securing  the  plant  to  the  floor  of  the  flume. 

The  instrumented  plant  in  the  small  flume  also  had  a  grommet,  but  was 
attached  to  a  smooth  aluminum  plate  rather  than  the  plywood  floor.  The  plate 
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was  15.2-cm  (6-in)  wide  by  30.5-cm  (12-in)  long  and  2.54-cm  (1-in)  thick.  The 
plate  provided  a  platform  by  which  to  measure  the  drag  force  produced  on  the 
plant.  A  hole  was  drilled  in  the  plate  and  a  shorter  grommet  was  used  because 
the  plate  was  not  as  thick  as  the  false  deck.  The  plant  was  inserted  through  the 
washer  and  the  grommet  screwed  to  the  plate  in  the  same  method  as  with  the 
other  plants. 

To  allow  placement  of  the  plate  into  the  flume  floor,  a  16.5-cm  (6-V2  in.)  by 
31.8  cm  (1244  in.)  rectangle  was  cut  in  the  center  of  the  floor  along  the  center 
line  of  the  flume.  Since  the  floor  was  3.81 -cm  ( 1  -Vi  in.)  thick,  1.27-cm-  04-in.-) 
diam  ball  bearings  were  placed  directly  on  the  smooth  steel  floor  where  the 
plywood  was  removed.  This  allowed  the  plate  to  move  smoothly  on  the  steel 
deck.  It  also  raised  the  top  of  the  plate  to  a  height  of  3.81  cm  (1 44  in.),  exactly 
flush  with  the  rest  of  the  floor.  This  prevented  the  water  from  striking  the  face  of 
the  plate  and  adding  to  the  measured  drag  force. 

The  same  strain  gage  used  in  the  large  flume  experiments  was  used  in  the 
small  flume.  It  was  placed  and  centered  directly  behind  the  aluminum  plate  to 
measure  the  drag  force  as  compression.  The  gage  and  connections  were  sealed  in 
waterproof  bags.  The  strain  gage  was  temperature  compensating  and  always 
zeroed  in  place  and  underwater.  The  calibration  of  the  gage  was  checked  before 
each  series  of  experiments. 

Elastic  bands  or  springs  were  attached  to  both  the  plate  and  the  plywood 
floor  immediately  downstream  and  to  the  sides  of  the  plate.  This  held  the  plate 
firmly  in  contact  with  the  strain  gage  and  centered  in  die  floor  cavity. 

Velocity  measurements  were  made  using  a  propeller  type  Ott  velocity  meter. 
Velocity  measurements  were  taken  at  the  center  of  the  leaf  mass  just  upstream  of 
the  instrumented  plant. 


Small  Sectional  Flume  Operating  Procedures 

Measurements  were  made  of  plant  dimensions  and  plant  characteristics 
before  each  series  of  experiments.  Plant  height  and  width,  leaf  size,  and  stem 
height  were  measured.  The  number  of  branches,  stems,  and  leaves  were  counted. 
The  diameter  of  stems  and  branches  was  recorded,  and  bending  characteristics 
were  measured. 

The  roots  of  the  plants  were  then  removed  and  attached  to  either  the  plywood 
floor  or  the  aluminum  plate.  Stop  logs  were  placed  to  a  height  of  0.914-m  (3  ft) 
at  the  downstream  end  of  the  flume.  This  allowed  the  flume  to  be  completely 
filled  and  the  strain  gage  set  to  0  to  compensate  for  any  buoyancy  effects. 

Each  plant  was  subjected  to  a  series  of  10  runs.  Each  run  was  at  an 
increasing  velocity,  ranging  from  approximately  0.076  to  2.43  m/sec  (0.25  to 
8  ft/sec.)  During  each  ran,  the  velocity  directly  upstream  of  the  plant  and  the 
compression  on  the  strain  gage  were  recorded.  This  velocity  was  taken  at  the 
center  line  of  the  effective  leaf  area.  As  velocity  increased,  the  velocity  probe 
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was  lowered  to  compensate  for  plant  bending.  This  insured  that  the  velocity  of 
each  run  was  being  recorded  at  the  vertical  center  line  of  the  leaf  mass.  The 
angle  the  plant  deflected  was  determined  from  marks  drawn  on  the  sidewalls  of 
the  flume.  Videotapes  were  taken  to  allow  for  more  detailed  observation  of  the 
plants  at  a  later  time. 

After  the  plant  was  subjected  to  10  different  velocities,  all  leaves  were 
removed.  The  plant  was  then  immediately  subjected  to  10  more  runs.  Velocity, 
drag,  and  deflection  data  were  again  recorded  as  previously  described. 
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3  Results  and  Analysis 


Resistance  Coefficients  from  Head  Loss 
Measurements 

Twenty-one  different  series  of  experiments  were  completed  in  the  large 
flume  using  different  plant  types,  plant  combinations,  plant  heights,  plant 
spacings,  flow  velocities,  and  depths.  Plant  characteristics  are  listed  in  Tables  1 
and  2  in  SI  and  non-SI  units,  respectively.  The  runs  were  videotaped  as  well  as 
photographed.  One  run  was  made  to  determine  the  bed  roughness  of  the  flume 
without  plants  on  each  bed  type  (clay  and  mortar). 

Tables  3  and  4  present  results  from  the  large-flume  homogeneous-plant¬ 
grouping  experiments  in  SI  and  non-SI  units,  respectively.  Tables  5  and  6  present 
the  results  from  large-flume  mixed-plant-grouping  experiments  in  SI  and  non-SI 
units,  respectively.  Water-surface  elevation,  average  depth,  and  discharge  were 
measured.  Average  velocity  was  determined  from  the  continuity  equation.  The 
average  roughness  coefficient  came  from  an  iterative  solution  of  the  backwater 
equation  in  which  calculated  water-surface  elevations  were  matched  to  measured 
water-surface  elevations.  In  the  backwater  equation,  Equation  2  was  used  to 
determine  n  and  Equation  3  to  determine  V*/V.  Average  energy  slope  is 
presented  in  Tables  3-6.  Flume  wall  effects  were  accounted  for  by  the  method 
advanced  by  Vanoni  and  Brooks  (1957). 

A  typical  velocity  profile  measurement  is  shown  in  Figure  5.  The  profile 
demonstrates  the  effect  of  leaf  mass  on  velocity.  The  plant  approach  velocity  is 
the  velocity  that  occurred  upstream  at  the  center  line  of  the  leaf  mass  of  the  plant. 
The  velocity  significantly  increased  below  the  leaf  mass.  Velocity  profile 
measurements  are  reported  in  Rahmeyer  et  al.  1996. 

Average  channel  velocities  between  0.914  and  1.22  m/s  (3  and  4  ft/sec)  were 
necessary  to  cause  either  the  leaves  to  break  off  the  plants  or  for  the  stems  to 
break.  These  velocities  also  caused  significant  movement  of  bed  material  in  the 
Phase  I  experiments.  It  is  possible  that  many  of  the  leaf  and  stem  failures  may 
have  been  due  to  the  impact  from  large  bed  material,  i.e.,  gravel  size  particles, 
being  transported  by  the  flow. 

Another  observation  was  that  shrubs  with  open  areas  beneath  the  primary 
leaf  mass  were  diverting  significant  amounts  of  flow  beneath  the  leaf  mass.  In 
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Figure  5.  Example  velocity  profile  for  an  experimental  run  with  dogwoods  (To 
convert  feet  to  meters,  multiply  by  0.3048) 


some  cases  velocities  below  the  leaf  mass  approached  surface  velocities  and  were 
sufficient  to  transport  even  the  largest  size  gravel  particles. 

The  Euonymus  plants  were  ground  cover  plants  with  leaves  extending  to  the 
bed.  These  plants,  when  used  in  a  typical  spacing,  left  areas  of  the  bed  exposed 
to  flow.  Measurable  scour  was  noted  in  these  areas  for  all  Phase  I  experiments. 
The  Euonymus  plant  experiments  were  stopped  when  scour  began  to  wash  the 
plants  away.  Only  the  wire  attached  to  the  plant  stems  kept  the  plants  from  being 
washed  downstream.  It  was  observed  that  local  scour  was  occurring  from  three- 
dimensional  vortices  that  appeared  to  be  similar  to  those  typically  associated  with 
scour  around  bridge  piers. 

Figures  6-11  demonstrate  the  effect  of  velocity  on  plant  deformation, 
sediment  transport,  and  scour. 


Calculation  of  Roughness  Coefficients 

The  hydraulic  roughness  and  the  Manning’s  coefficient  n  for  plant  resistance 
were  calculated  by  using  an  initial  estimate  of  a  total  Manning's  roughness 
coefficient  and  then  adjusting  the  n  value  to  best  fit  the  gradually  varied 
backwater  curve  of  measured  water-surface  elevations  along  the  experimental 
section.  Equation  8  was  the  equation  used  to  fit  the  backwater  curve. 
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Figure  7.  Plants  at  low  flow 
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Figure  8.  Plants  at  moderate  flow 


Figure  9.  Plants  with  sediment  transport 
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dy_ 

dx 


(8) 


rs0-s,  'i 

Here  dy/dx  is  the  unit  change  in  slope  of  the  water  surface;  S0  is  the  slope  of  the 
bed;  Sf  is  the  slope  of  the  energy  line;  and  Fris  the  Froude  number.  Sf  is 
calculated  from  the  Manning's  equation  (Equation  1)  using  the  estimate  of 
Manning’s  n,  the  mean  velocity,  V,  calculated  from  the  continuity  equation,  and 
the  hydraulic  radius,  Rh.  The  Froude  number  was  calculated  from  Equation  9. 


V 


(9) 


The  Manning’s  n  for  the  vegetation  and  the  bed  was  then  iteratively  solved 
using  a  trial  and  error  process  until  the  shape  of  the  backwater  curve  predicted  by 
Equation  8  was  the  same  as  the  measured  curve  of  the  actual  water  surface. 
Figure  12  is  an  example  of  the  backwater  curve  fit  for  a  run  with  a  total 
Manning’s  n  of  0.062. 


Figure  12.  Typical  backwater  curve  for  experimental  runs  (To  convert  feet  to 
meters,  multiply  by  0.3048) 


From  the  total  Manning’s  n  the  value  of  the  bed  and  plant  roughness,  was 
determined.  The  first  step  in  this  determination  was  to  convert  the  total  n  to  a 
Darcy-Weisbach  friction  factor,/,  by  Equation  10. 
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The  coefficient  of  friction  for  the  bed  and  plants,  /*,  was  determined  using  a 
correction  for  flume  wall  effects.  The  coefficient  of  friction  for  the  walls,/*,,  was 
determined  from  Equation  11.  This  equation  was  regressed  for  this  study  to  fit 
the  correction  figure  presented  by  Vanoni  and  Brooks  (1957). 
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0.274367 
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Re  is  the  Reynolds  number.  Equation  1 1  was  a  power  fit  regression  with  an  R2 
of  99.98  percent.  The  friction  factor  for  the  bed  and  plants,/*,  was  then 
calculated  with  Equation  12. 


ft  =  /  -  fw) 

b 


(12) 


Here  b  is  the  width  of  the  channel,  and  Y0  is  the  flow  depth.  The  hydraulic  radius 
associated  with  the  bed  and  plants,  Rb,  was  determined  by  Equation  13. 


Rb_  _  Rw_  _ 

fb  fw  f 


(13) 


Rw  is  the  hydraulic  radius  for  the  walls;  and  Rh  is  the  total  hydraulic  radius. 
Equations  12  and  13  are  from  Vanoni  and  Brooks  (1957).  Finally,  the  Manning’s 
coefficient  n*  for  the  bed  and  vegetation  roughness  was  calculated  using  Rb  in 
Equation  1. 


The  coefficient  n*  is  the  resistance  of  both  the  bed  and  vegetation  roughness. 
Equation  14  can  be  used  to  calculate  the  resistance  coefficient  nveg  for  the  net 
resistance  of  the  vegetation. 


n  veg  n  b  n  0 


(14) 


where  nveg  is  the  Manning’s  coefficient  for  vegetation;  n*is  the  bed  and  vegeta¬ 
tion  resistance;  and  n<>  is  the  bed  roughness.  The  value  for  n^  for  both  the  clay 
and  mortar  beds  (corrected  for  wall  effects)  was  determined  to  be  approximately 
0.020.  V./V  was  found  to  have  a  value  of  0.069.  Resistance  coefficients 
reported  in  this  report  show  the  combined  value  for  the  bed  and  the  plants.  This 
combined  value  is  typically  reported  in  field  investigations.  Roughness 
coefficients  reported  herein  may  be  reduced  to  strictly  vegetation  roughness 
coefficients  by  subtracting  0.020  from  reported  Manning’s  n  values  and  0.069 
from  reported  V./V  values. 
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Measured  Drag  Forces 

Drag  forces  were  measured  in  both  the  large  and  small  flumes.  More  than 
100  experiments  were  conducted  in  the  small  flume  using  10  different  plant 
types,  and  more  than  20  experiments  were  conducted  in  the  large  flume  using 
four  plant  types.  Measured  approach  velocities  and  drag  forces  are  shown  in 
Tables  7  and  8  for  the  small  and  large  flumes,  respectively.  Data  from  four 
different  dogwood  plants  are  plotted  in  Figure  13  showing  the  repeatability  of 
drag  force  measurements  between  the  large  and  small  flumes.  This  is  important 
because  it  demonstrated  that  the  experimental  data  from  the  small  flume  could  be 
directly  compared  to  the  plants  and  resistance  coefficients  determined  in  the  large 
flume.  In  Figure  13,  note  that  the  drag  force  increases  linearly  with  velocity 
instead  of  with  the  square  of  velocity  as  one  would  expect  from  the  drag  force 
equation  (15).  This  occurred  in  the  experiment  because  the  drag  coefficient  and 
blockage  area  were  reduced  as  the  plant  was  streamlined  by  the  increasing 
velocities. 


Fd  = 


pcdav; 


(15) 


In  Equation  15,  FD is  the  drag  force,  p  is  the  density  of  water,  CD  is  the  drag 
coefficient,  Vpis  the  approach  velocity,  and  A,  is  the  blockage  area  of  an 
individual  plant. 


Figure  1 3.  Plant  approach  velocity  versus  drag  force  (T o  convert  feet  to  meters, 
multiply  by  0.3048)  [To  convert  pounds  (force)  to  newtons,  multiply 
by  4.448222] 
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The  small  flume  had  a  large  plastic  window  through  which  plant  distortion 
could  be  viewed  and  measured.  It  was  observed  that  the  plants  easily  bent  with 
the  flow,  and  the  leaf  mass  trailed  downstream  forming  a  streamlined,  almost 
teardrop-shaped  profile.  The  leaf  mass  changed  with  velocity  and  became  more 
streamlined  with  increased  velocity.  This  observation  explains  the  significant 
decrease  in  resistance  with  increase  in  velocity.  The  changing  shape  of  the  leaf 
mass  means  that  the  roughness  coefficient  will  change  with  velocity  and  that  the 
assignment  of  a  constant  roughness  coefficient  to  determine  a  stage-discharge 
curve  would  be  invalid  and  produce  significantly  incorrect  results. 

If  plant  resistance  were  constant  with  increasing  velocities,  a  plot  of  velocity 
versus  drag  force  would  appear  as  a  smooth  exponentially  increasing  curve.  A 
typical  curve  from  the  data  is  shown  in  Figure  14.  In  this  figure  the  drag  varies 
almost  linearly  with  velocity  as  the  leaf  mass  continues  to  streamline.  At  a 
velocity  of  about  1.2  m/sec  (4  ft/sec)  the  leaf  mass  has  reached  its  streamlining 
limit  and  the  curve  begins  to  take  on  the  expected  exponential  form.  Above  this 
limiting  velocity  it  would  be  appropriate  to  assign  a  constant  roughness 
coefficient. 


Figure  14.  Drag  force  versus  velocity  (To  convert  feet  to  meters,  multiply  by 

0.3048)  [To  convert  pounds  (force)  to  newtons,  multiply  by  4.448222] 


Development  of  Resistance  Methodology 

The  methodology  developed  in  this  report  was  based  on  the  premise  that  the 
flow  resistance  due  to  vegetation  on  a  floodplain  is  equal  to  the  sum  of  the  total 
drag  forces  produced  by  that  vegetation.  Kadlec  (1990)  presented  such  a 
hypothesis  and,  assuming  that  the  drag  forces  on  the  bed  are  negligible,  gave  the 
following  conceptual  equation  relating  shear  stress  to  drag  force. 


_  y  Fp 

area 


FdM 


(16) 
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where 


area  =  total  bed  surface  area  of  interest 


T0  =  shear  stress  (  t0  =  y  Rh  S  ) 


M  =  plant  density 

The  shear  velocity  V*  is  related  to  shear  stress,  and  a  commonly  used 
resistance  coefficient  that  is  associated  with  shear  stress  is  V»/V  (Equation  2). 


V*_  _ 

V  V 


l PV 2 


(2) 


The  shear  velocity  and  shear  stress  can  then  be  related  to  drag  force  by  using 
Equation  17: 


V* 


(17) 


Using  Equations  15  through  17,  an  equation  relating  the  drag  coefficient,  CDito 
the  resistance  can  be  developed. 


Yi 
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( At  MCd  \ 


0.5 


(18) 


The  blockage  area,  Ai;  of  a  plant  with  dense  foliage  is  approximated  by 
multiplying  the  effective  plant  height  times  the  effective  plant  width  (H'  x  We). 
This  effective  area  is  the  equivalent  rectangular  area  of  the  leaf  mass  discounting 
small  stems  that  are  not  part  of  the  average  leaf  mass.  For  plants  that  have  voids 
in  their  leaf  mass  or  few  leaves  with  a  large  number  of  stems,  A*  is  the  rectangu¬ 
lar  area  equivalent  to  the  net  frontal  blockage.  For  example,  plants  without 
leaves  would  have  a  blockage  area  equal  to  the  number  of  stems  times  the  stem 
diameter  times  the  stem  length.  Blockage  areas  for  the  laboratory  experiments 
were  determined  from  digital  photographs  of  the  plant  against  a  white  back¬ 
ground  marked  with  grid  lines.  H'  is  the  actual  height  of  the  undistorted  leaf 
mass  and  We  is  the  effective  width  that  produces  the  measured  Ai. 

It  has  been  established  that  the  drag  coefficient  for  a  rigid  body  is  not  a 
constant  and  varies  with  Reynolds  number,  Re.  The  Reynolds  number  used  in 
this  study  is  based  on  the  length  variable  of  hydraulic  radius,  Rh,  and  the  mean 
channel  velocity.  For  flexible  plants,  the  drag  coefficient  is  a  function  of  a 
number  of  factors  as  shown  in  Equation  19. 

CD  =f  (Re,  Ya,  H,  plant  type,  plant  shape,  plant  flexibility,  M)  (19) 
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The  experimental  data  from  both  the  large  and  small  flumes  were  used  to 
determine  appropriate  dimensionless  parameters  to  define  the  drag  force.  These 
experiments  were  conducted  for  a  large  matrix  of  variables  including  Y0,  V,  plant 
type,  leaf  density,  plant  density,  plant  shape,  plant  size,  and  blockage  area.  The 
runs  were  made  in  a  sequence  so  that  each  variable  could  be  evaluated  by  keep¬ 
ing  the  other  variables  constant.  It  was  found  that  the  drag  coefficient  decreased 
with  an  increase  in  velocity,  depth,  plant  density,  plant  flexibility,  and  plant  spac¬ 
ing.  Drag  coefficient  or  resistance  could  not  be  related  solely  to  flow  conditions, 
leaf  density,  or  plant  blockage  because  of  the  flexibility  of  the  plants.  For 
example,  different  plants  with  the  same  size  leaves  and  blockage  had 
significantly  different  resistance  depending  upon  how  much  the  plants  deformed 
and  bent  with  flow. 

Dimensional  analysis  was  used  to  aid  in  the  selection  of  dimensionless 
parameters  that  could  relate  drag  or  resistance  to  flow  and  plant  variables.  The 
four  parameters  that  were  found  to  have  a  significant  effect  were  as  follows: 

a.  Ratio  of  the  flow  drag  force  to  the  forces  resisting  plant  distortion. 

b.  Ratio  of  the  flow  depth  to  the  plant  height. 

c.  Blockage  of  the  plants  to  the  flow  on  the  channel  bottom. 

d.  Reynolds  number. 

The  last  three  parameters  are  corrections  or  modifiers  to  the  ratio  parameter  of 
drag  force  to  the  force  resisting  plant  deformation.  This  ratio  parameter  also 
incorporates  plant  stiffness  or  flexibility.  These  parameters  are  shown  in 
Equation  20. 


CD  or  — 
D  V 


=  f 


6V1  A-  Y 

- - 2-,  -2-,  M  A  ,  Re 

EA  H  ^  e 


(20) 


When  there  are  several  stems  emerging  from  the  base  of  the  plant,  the 
stiffness  modulus,  Es,  is  determined  for  a  single  stem.  The  stem  area,  As,  is  the 
sum  of  all  the  stem  areas.  The  stem  area  for  a  plant  with  multiple  stems  is  thus 
calculated  as  the  number  of  stems  times  tc  times  the  stem  diameter  squared 
divided  by  4. 


Resistance  Equation  for  Submerged  Vegetation 

The  results  from  the  large  flume  experiments  were  analyzed  to  determine  the 
regression  of  the  variables  of  Equation  20  for  submerged  vegetation.  The 
regression  analysis  found  that  log  and  polynomial  relationships  gave  a  poor  data 
fit  while  a  power  relationship  had  very  good  results.  Equations  21  and  22  were 
found  to  fit  the  data  with  a  regression  coefficient  of  R2  =  96  percent  and  a 
maximum  scatter  of  15  percent  for  predicted  values  of  V./V  with  measured 
values.  The  parameters  in  the  equations  were  modified  to  allow  a  direct  solution 
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for  resistance  (for  a  given  depth)  by  combining  the  original  parameters  with 
Manning’s  equation  and  the  equation  for  shear  velocity.  This  modification  and 
combination  of  equations  resulted  in  Equations  21  for  shear  velocity  and 
Equation  22  for  Manning’s  n. 
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It  is  important  to  note  that  the  plant  characteristics  H,  Aj,  and  As  are  the 
initial  characteristics  of  the  plants  without  the  effects  of  flow  distortion.  During 
the  experiments,  it  was  observed  that  since  the  plants  bent  with  flow, 
submergence  occurred  at  flow  depths  less  than  80  percent  of  the  plant  height. 
Equations  21  and  22  are  to  be  applied  only  for  submerged  flow  defined  by  Y0 
>0.8  H.  Equations  23  and  24  are  for  partially  submerged  flow  with  Y0  <0.8  H. 
Both  sets  of  equations  converged  to  approximately  the  same  result  at  the  flow 
depth  Y0  =0.8  H. 


Resistance  Equation  for  Partially  Submerged 
Vegetation 

The  data  for  partially  submerged  vegetation  were  analyzed  to  determine  the 
regression  of  the  variables  of  Equation  20.  The  regression  analysis  again  found 
that  a  log  relationship  gave  a  poor  fit  of  data  while  a  power  relationship  produced 
very  good  results.  Equations  23  and  24  were  found  to  fit  the  data  with  a 
regression  coefficient  of  R2=  85  percent  and  a  maximum  scatter  of  18  percent  for 
predicted  values  of  V*/V  compared  to  measured  values.  These  equations  again 
allow  direct  solution  for  resistance  if  the  flow  depth  is  known. 
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The  blockage  area  in  Equations  23  and  24  was  changed  to  an  effective  area, 
A since  only  a  portion  of  die  leaf  mass  produces  blockage  under  partially 
submerged  flow  conditions.  The  effective  blockage  area  can  be  approximated 
using  Equation  25  if  the  actual  geometry  of  the  plant  and  leaf  mass  has  not  been 
measured  (see  Figure  2). 
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A*  =  [Y0-{H-H')]We 


(25) 


The  analysis  of  data  and  the  regression  fit  of  Equations  21  through  24 
included  many  other  parameters  and  ratios.  The  equations,  parameters,  and 
methods  developed  by  other  researchers  for  a  combined  density  and  blockage  of 
heavy  ground  cover  and  grasses  did  not  produce  satisfactory  results.  These 
methods  included  those  developed  by  Kowen  and  Li  (1980),  Ree  and  Crow 
(1977),  and  other  methods.  It  should  be  noted  that  some  of  the  methods 
evaluated  were  not  developed  with  shrubs  in  mind.  In  the  case  of  Ree  and  Crow, 
agricultural  crops  were  the  focus  of  the  methodology.  Many  of  the  equations 
developed  by  other  researchers  were  evaluated,  but  none  proved  to  be 
satisfactory  in  the  prediction  of  roughness  values  for  shrubs.  The  results  of  this 
study  emphasize  that  the  plant  stiffness  modulus,  Es,  must  be  considered  to 
obtain  a  satisfactory  prediction  of  roughness.  The  definition  and  method  to 
determine  Es  are  discussed  in  a  following  section. 

Equations  21  through  24  also  include  plants  with  multiple  stems.  The 
blockage  area  A;  is  for  an  individual  or  average  plant,  the  plant  density  is  the 
number  of  plants  (not  stems)  per  unit  area,  and  As  is  the  sum  of  the  cross- 
sectional  area  of  all  of  the  stems  of  an  individual  average  plant.  Figure  15  shows 
the  correlation  of  calculated  V./V  with  the  observed  data  from  the  flumes  for 
submerged,  partially  submerged,  multiple  plant  species. 
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Multiple  Plant  Combinations 

Six  combinations  of  plant  types  were  evaluated.  The  combinations  of  species 
were  selected  to  represent  typical  plant  groups  found  in  different  ecosystems. 
Equations  21  through  24  worked  equally  well  for  plant  combinations  when 
average  plant  characteristics  were  used.  The  average  plant  characteristics  were 
obtained  by  weighting  individual  plant  characteristics  by  the  number  of  plants 
per  unit  area  or  densities  of  each  type  of  plant.  The  comparisons  of  calculated 
and  observed  data  for  multiple  plant  groupings  are  shown  in  Figure  15. 

The  purpose  of  the  weighted  averages  is  to  formulate  the  average  shear  stress 
created  by  the  plant  combinations.  Each  plant  group  then  will  have  an  average 
blockage  area  or  effective  blockage  area,  an  average  modulus  of  plant  stiffness, 
an  average  total  plant  stem  area,  an  average  plant  height,  and  an  average  effective 
plant  height.  A  weighted  average  for  the  plant  groups  is  then  based  on  the  ratio 
of  the  plant  density  of  each  plant  type  divided  by  the  total  plant  density  of  all  the 
plants  and  plant  types.  Equations  21  through  24  do  not  use  an  average  plant 
density,  but  use  the  total  or  combined  density  of  all  of  the  plants.  The  method  for 
combining  these  densities  is  shown  in  the  example  that  follows. 


The  plant  characteristics  are  determined  by  weighting  the  individual  plant 
characteristics  according  to  their  relative  density  in  the  area  of  interest.  The 
weighted  values  are  then  summed  to  obtain  the  combined  plant  characteristic  as 
shown  in  Equations  26-32. 
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Using  Equation  29,  the  average  plant  height  for  an  area  with  a  group  of  three 
plants  with  heights,  H,  of  0.5, 1.0,  and  2.0  ft;  and  densities,  M,  of  0.25, 0.50,  and 
0.20  plants  /  ft2;  would  be  determined  as  follows: 


Havg  =  Hi  x  Mi  /  M total  +  H2  x  M2  /  +  H3  x  M3  /  Mtotai 


Havg  =  0.5  x  0.25  /  (0.25+0.5+0.2)  +1.0  x  0.5  /  0.95  +  2.0  x  0.2  /  0.95 
Havg  =  0.132  +  0.526  +  0.421  =  1.079  ft 


Stiffness  Modulus 


The  modulus  of  plant  stiffness,  Es,  is  critical  to  the  calculation  of  resistance 
because  of  the  flexibility  of  the  plants  and  the  deformation  of  leaf  masses  due  to 
the  flow  forces.  The  modulus  of  plant  stiffness  is  calculated  by  Equation  33. 
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The  data  necessary  to  use  Equation  33  is  obtained  by  measuring  the  force, 
F45,  necessary  to  bend  the  plant  to  an  angle  of  45  deg.  The  45-deg  angle  is 
measured  from  the  initial  vertical  position  to  the  stem  or  leaf  mass  at  the  point 
where  the  force  is  measured,  i.e.,  at  H/2  as  shown  in  Figure  16. 

I  is  the  second  area  moment  of  inertia  calculated  for  a  circular  shape  (I  = 
7tDs4/64).  The  stem  diameter  Ds  is  measured  at  a  height  of  H/4  above  the  ground. 

Data  were  collected  both  in  the  laboratory  and  in  the  field  to  determine  a 
relationship  that  defined  plant  stiffness.  Freeman  (1997)  collected  data  for  five 
types  of  willows  in  floodplains  and  on  sand  bars  to  determine  if  stiffness  in  the 
field  could  be  predicted  from  plant  size  parameters  such  as  stem  diameter  and 
plant  height  to  reduce  the  number  of  parameters  that  must  be  collected  to 
determine  the  plant  stiffness  modulus.  Data  collected  included  samples  from 
Salix  exigua  willows  in  Utah  and  Idaho,  Salix  lasiandra,  Salix  lamonii,  a  wild 
rose  bush  common  to  the  area,  and  young  cottonwood  trees  growing  on  sandbars. 
He  also  noted  in  his  data  collection  efforts  that  plant  stiffness  was  measurably 
different  in  the  upstream  and  downstream  directions  in  streams  subject  to  long 
periods  of  high  water  (i.e.,  snow  melt).  Where  the  plants  were  not  subject  to 
velocities  high  enough  to  keep  the  plant  bent  and  /or  deformed  for  prolonged 
periods  of  time  this  difference  in  the  stiffness  modulus  did  not  seem  to  exist  or 
was  not  noticeable. 
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Figure  16.  Methodology  for  measurement  of  plant  stiffness  for  calculation  of  Es 
in  the  field  for  plants  with  effective  height  of  leaf  mass  approximately 
equal  to  the  plant  height 

The  research  performed  in  the  laboratory  and  in  the  field  indicated  that 
the  stiffness  modulus  can  be  estimated  from  the  relationship  of  Es  to  the  ratio  of 
H/Ds.  The  analysis  of  measurements  made  in  the  field  and  in  the  laboratory  led 
to  the  development  of  Equations  34  and  35  to  explain  the  relationship  between 
H/Ds  and  Es.  The  relationship  between  the  data  observed  in  Freeman’s  field 
measurements  (Freeman  1997,  Freeman,  et  al.  1998)  and  the  values  predicted  by 
Equation  34  is  shown  in  Figure  17.  Equation  34  (shown  in  Figure  17  as 
“Rahmeyer  Predicted”)  was  developed  based  on  laboratory  observations  and 
gives  the  modulus  in  pounds  per  square  foot  while  Equation  35  gives  the  value  in 
newtons  per  square  meter.  It  must  be  cautioned  that  the  fit  of  Equations  34  and 
35  have  a  regression  (R2)  of  less  than  90  percent,  and  the  scatter  is  significant  as 
shown  in  Figure  17. 
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Actual  field  measurements  of  Es  are  recommended  where  possible.  The  stiffness 
modulus  can  also  be  estimated  from  measured  values  of  similar  plants.  Since  the 
stiffness  modulus  varies  depending  on  the  plant  size,  it  was  determined  that  if  the 
calculated  modulus  for  a  particular  plant  size  was  divided  by  (H/Ds)1 5,  the 
stiffness  modulus  became  independent  of  plant  size  and  one  value  could  be  used 
for  all  plant  sizes.  Thus,  to  calculate  the  plant  stiffness  modulus  for  Alder,  the 
value  from  Table  9  is  multiplied  by  (H/Ds)1 5  which  gives  the  stiffness  modulus 
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Figure  17.  Measured  versus  calculated  plant  stiffness  modulus,  Es 

for  a  particular  size  plant  in  pounds  per  square  foot.  The  exponent  for  the  term 
H/Ds  was  determined  to  remove  most  effects  of  plant  size  from  the  Plant 
Stiffness  Modulus. 
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A  total  of  20  different  plant  species  were  evaluated  in  either  a  large  2.44-m- 
(8-ft-)  wide  flume  or  a  small  0.46-m-  (1.5-ft-)  wide  flume  to  determine  flow 
resistance  and  drag  force.  More  than  220  experiments  were  conducted.  Fifteen 
homogeneous  plant  groupings  were  evaluated  in  the  large  flume.  Six  multiple 
plant  groupings  were  evaluated  in  the  large  flume.  Vegetation  was  evaluated 
under  both  submerged  and  partially  submerged  conditions.  Velocity,  depths, 
plant  density,  plant  dimensions,  and  plant  types  were  varied  in  the  experiments. 
The  range  of  experimental  conditions  used  to  develop  the  regression  equations 
were  as  follows: 

a.  Flow  depths  from  0.4  to  1.4  m  (1.3  to  4.7  ft). 

b.  Average  flow  velocities  from  0.15  to  1.1  m/s  (0.5  to  3.6  ft/s). 

c.  Measured  resistance  V*/V  from  0.13  to  0.43  and  n  from  0.04  to  0.14. 

d.  Plant  heights  from  0.20  to  1.52  m  (0.66  to  5  ft). 

e.  Plant  widths  from  0.076  to  0.91  m  (0.25  to  3  ft). 

f  Plant  densities  from  0.53  to  13  plants  /  m2  (0.05  to  1.2  plants  /  ft2). 

g.  Plant  modulus  of  stiffness  from  5.3  x  107  to  4.8  x  109  N/m2  (1.1  x  106to 
1.0  x  lO^bfrft2). 

h.  Reynolds  numbers  from  1.4  x  105  to  1.6  x  106. 

An  important  observation  made  during  the  flume  studies  was  that  the  plant 
leaf  mass  trailed  downstream  forming  a  streamlined,  almost  teardrop-shaped 
profile.  The  leaf  mass  shape  changed  with  velocity  and  became  more 
streamlined  with  increasing  velocity.  The  effect  of  this  phenomenon  was  a 
significant  decrease  in  the  drag  coefficient  and  resistance  coefficient  with 
velocity.  On  the  other  hand,  resistance  increased  with  depth  for  partially 
submerged  plants  as  the  blockage  area  increased  with  depth  until  the  plants  were 
submerged.  The  transition  between  submerged  and  partially  submerged  flow 
occurred  at  a  depth  of  about  80  percent  of  the  undeflected  plant  height. 
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Another  observation  made  during  the  study  was  that  the  leaf  mass  or  foliage 
canopy  diverted  flow  beneath  the  canopy.  The  bottom  flow  resulted  in 
significant  velocities  along  the  channel  bed  causing  general  scour  and  increased 
sediment  transport.  The  bed  velocities  were  sufficient  to  transport  and  move  the 
largest  sizes  of  gravel  found  in  the  flume  bed. 

The  hydraulic  roughness  of  a  vegetated  channel  was  shown  to  be  a  function 
of  the  stiffness  of  the  plants  growing  in  the  channel,  the  depth,  velocity,  and 
hydraulic  radius  of  the  channel,  plant  density,  and  frontal  area  of  the  plant 
obstructing  the  flow.  It  was  determined  that  the  roughness  can  be  calculated 
directly  if  the  depth  of  flow  is  known.  The  roughness  can  be  determined  in  terms 
of  Manning’s  n,  Chezy  C,  or  the  shear  velocity  ratio,  V»/V. 

Regression  equations  were  developed  for  submerged  vegetation  and  found  to 
fit  the  data  with  a  regression  coefficient  of  R2  =  96  percent  and  a  maximum 
scatter  of  15  percent  for  predicted  values  of  V*/V  with  measured  values.  The 
parameters  in  the  equations  were  modified  to  allow  a  direct  solution  for 
resistance  (for  a  given  depth)  by  combining  the  original  parameters  with 
Manning’s  equation  and  the  equation  for  shear  velocity.  This  modification  and 
combination  of  equations  resulted  in  Equations  21  for  shear  velocity  and 
Equation  22  for  Manning’s  n. 
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It  is  important  to  note  that  the  plant  characteristics  H,  A;,  and  As  are  the 
initial  characteristics  of  the  plants  without  the  effects  of  flow  distortion. 
Equations  21  and  22  are  to  be  applied  only  for  submerged  flow  defined  by 
Y0  >0.8  H. 
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The  experimental  data  for  partially  submerged  vegetation  were  analyzed  to 
develop  regression  Equations  23  and  24.  These  equations  were  found  to  fit  the 
data  with  a  regression  coefficient  of  R2=  85  percent  and  a  maximum  scatter  of 
18  percent  for  predicted  values  of  V*/V  compared  to  measured  values.  These 
equations  again  allow  direct  solution  for  resistance  if  the  flow  depth  is  known. 
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The  blockage  area  in  Equations  23  and  24  was  changed  to  an  effective  area.  A,*, 
since  only  a  portion  of  the  leaf  mass  produces  blockage  under  partially  sub¬ 
merged  flow  conditions.  The  effective  blockage  area  can  be  approximated  using 
Equation  25  if  the  actual  geometry  of  the  plant  and  leaf  mass  has  not  been 
measured. 

A,'  =  (25) 

The  resistance  coefficients  predicted  by  Equations  21  through  24  represent 
the  combined  resistance  of  the  bed  and  the  plants.  Resistance  coefficients  due 
only  to  vegetation  must  be  calculated  by  subtracting  the  bed  resistance.  In  these 
experiments  the  Manning’s  bed  resistance  coefficient  was  found  to  be  0.02  and 
V*/V  for  the  bed  was  found  to  be  0.069. 

The  modulus  of  plant  stiffness,  Es,  is  critical  to  the  calculation  of  resistance 
because  of  the  flexibility  of  the  plants  and  the  deformation  of  leaf  masses  due  to 
the  flow  forces.  The  research  performed  in  the  laboratory  and  in  the  field 
indicated  that  the  stiffness  modulus  can  be  estimated  from  the  relationship  of  Es 
to  the  ratio  of  H/Ds.  The  analysis  of  measurements  made  in  the  field  and  in  the 
laboratory  led  to  the  development  of  equations  to  explain  the  relationship 
between  H/Ds  and  Es.  The  equations  had  a  regression  (R2)  of  less  than 
90  percent.  Actual  field  measurements  of  Es  are  recommended  where  possible. 

The  stiffness  modulus  can  also  be  estimated  from  measured  values  of  similar 
plants.  Since  the  stiffness  modulus  varies  depending  on  the  plant  size,  it  was 
determined  that  if  the  calculated  modulus  for  a  particular  plant  size  was  divided 
by  (H/Ds)1'5,  the  stiffness  modulus  became  independent  of  plant  size  and  one 
value  could  be  used  for  all  plant  sizes.  Stiffness  modulus’  of  plants  used  in  this 
study  are  provided  in  Table  9. 
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Table  1 

Dimensions  and  Characteristics  of 


Plant 
Height  H, 
m 


Plants  (SI  Units) 


Effective 

Height, 

H, 


Blockage 
Area  A 


Stem 

Number 


Comus  Stolonifera  0.51 
Flaviramea 


0.71 


Purpleleaf 

Euonymus 


Red  Twig 
Doawood 


Service  Be 


Mulefat 


Alder 


Euonymus  Fortunei  0.20 
Colorata 


Comus  Sericea  0.97 


Amelanchier 


Salt  Cedar 


Black  Willow 


Red  Willow 


Purpleleaf 

Eunonvus 


Artie  Blue 
Williow 


Norway  Maple 


Common  Privet 


Blue  Elderberry 


French  Pink 
Pussywillow 


Sycamore 


Staghorn 

Sumac 


Baccharis  Glutmosa  0.97 


I  Alnuslncana  0.76 


Sambucus  Mexicana  0.97 


Tamarix  s. 


Salix  Niqra 


Salix  Monticola 


Comus  Stolonifera 
Flaviramea 


Euonymus  Fortunei  0.20 
Colorata 


Salix  Purpurea  Nana 


Acer  Platenoides 


nsmE 


0.71 


0.81 


Sambucus 

Canadensis 


Salix  Caprea 
Penduia 


Platenus  Acer  Ifolia  10.91 


Prunis  Besseyi 


Rhus  Typhina 


Salix  exigua 


Pacific  Willow  I  Salix  lasiandra 


Salix  Lemonii 


Rosa  spp. 


ItSE 


Small  Flume 


0.229  0.33 


0.33 

0.076  ! 

0.51 

0.181 

0.20 

0.052 

0.76 

0.368 

0.51 

0.090 

0.61 

0.155 

0.51 

0.039 

0.70 

0.107 

0.91 

0.697 

1.37 

0.836 

1.22 

0.372 

0.61 

0.093 

1.22 

1.115 

0.0063 


0.0095 


0.0268 


0.0316 


0.0189 


0.0095 


0.0254 


0.0063 


0.0126 


0.0126 


0.0126 


0.0252 


10.0101 


0.0084 


0.0126 


0.015 


0.0095 


12 


3.210 


0.526 


4.140 


10.200 


I  EE 


.600 


29.900 


15.950 


17.000 


13.100 


1.500 


14.500 


3.410 


19.100 


3.940 


27.500 


28.800 


Table  2 

Dimensions  and  Characteristics  of  Plants  (Non-Si 


Units) 


Plant 

Common 

Height  H, 

Name 

Scientific  Name 

ft 

Blockage 
Area  A 
Ft2 


Stem  Elasticity 

Number  Es,  lbf/fr 


Comus  Stolonifera  1 .67 
Flaviramea 


Sambucus  2.33 

Racemosa 


Purpleleaf  Euonymus  Fortune!  0.67 

Euonvmus  Colorata 


Comus  Sericea  3.18 


Service  Berry  I  Amelanchier  2.33 


Comus  Stolonifera  2.33 
Flaviramea 


Bacchans  Glutinosa  3.18 


Alnuslncana  2.50 


Sambucus  Mexicans  3.18 


Salt  Cedar  I T amarix  s 


Black  Willow  Salix  Niara 


Red  Willow 


Salix  Monticola 


Common 

Privet 


Comus  Stolonifera  1 .67 
Flaviramea 


Euonymus  Fortunei  0.67 
Colorata 


Salix  Purpurea  Nana  1 .84 


Acer  Platenoides  2.33 


Ligustrum  Vulgare  2.67 


Sambucus 

Canadensis 


French  Pink 
Pussywillow 


Sycamore  Platenus  Acer  Ifolia  1 3.00 


Prunis  Besseyi  2.43 


Rhus  Typhina 


Salix  exigua 


Staghorn 

Sumac 


Sand  Bar 
Willow 


Pacific  Willow  Salix  lasiandra 


Lemon's  Salix  lemonii 
Willow 


Wild  Rose  Rosa  spp. 
Bush 


Large  Flume 


0.750  1.08 


Small  Flume 


0.818 

0.0313 

0.560 

0.0208 

1.669 

1.001 

0.0417 

1.873 

0.0417 

1.997 

0.0833 

0.700 

0.0625 

1.831 

0.0333 

0.829 

0.0278 

0.829 

0.0417 

7.09 

0.0492 

21.31 

0.0558 

4.09 

0.0427 

0.0230 

6.706 


1.099 


8.648 


1 

6.706 

2 

8.648 

1 

2.486 

1 

39.900 

1 

8.231 

1 

0.549 

1 

2.319 

1 

57.448 

1 

60.163 

1 

10.612 

1 

180 

1 

207  j 

1 

180 

1 

272 

Table  3 
Summary  of 


Large  Flume  Results  with  Homogeneous  Groupings  (SI  Units) 


Bed 

Hydraulic 


Energy 

Slope  Average 


0.388 

0.00013 

0.562 

0.069 

0.0200 

0.209 

0.00002 

0.016 

0.884 

0.064 

0.0200 

0.591 

0.00015 

0.016 

1.011 

0.069 

0.0220 

0.366 

0.00053 

0.046 

1.202 

0.216 

0.0710 

0.610 

0.00124 

0.042 

1.184 

0.198 

0.0650 

0.750 

0.00184 

0.040 

1.059 

0.185 

0.0590 

0.482 

0.00119 

0.047 

0.902 

0.213 

0.0670 

0.588 

0.00140 

0.043 

0.971 

0.196 

0.0620 

0.689 

0.00163 

0.040 

0.991 

0.183 

0.0580 

0.878 

0.00582 

0.048 

0.521 

0.197 

0.0560 

0.991 

0.00477 

0.041 

0.688 

0.181 

0.0540 

1.091 

0.00418 

0.038 

0.843 

0.170 

0.0530 

0.765 

0.00102 

0.031 

1.232 

0.145 

0.0480 

0.924 

0.00165 

0.031 

1.056 

0.142 

0.0460 

1.058 

0.00693 

0.040 

0.499 

0.174 

0.0500 

0.750 

0.00496 

0.042 

0.421 

0.191 

0.0530 

0.294 

0.00030 

0.042 

1.134 

0.195 

0.0640 

0.479 

0.00063 

0.035 

0.918 

0.157 

0.0500 

0.589 

0.00085 

0.034 

0.989 

0.154 

0.0490 

0.304 

0.00043 

0.045 

0.908 

0.204 

0.0640 

0.518 

0.00125 

0.040 

0.676 

0.176 

0.0530 

0.614 

0.00110 

0.033 

0.735 

0.145 

0.0440 

0.00123 

0.032 

0.793 

0.141 

0.0430 

0.769 

0.00167 

0.033 

0.767 

0.146 

0.0450 

0.862 

0.00199 

0.031 

0.702 

0.136 

0.0410 

0.945 

0.00191 

0.030 

0.849 

0.133 

0.0410 

0.319 

0.00041 

0.045 

1.120 

0.209 

0.0680 

0.420 

0.00055 

0.040 

1.122 

0.186 

0.0600 

0.669 

0.00159 

0.042 

1.063 

0.195 

0.0630 
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Table  4 

Summary  of  Large  Flume  Results  with  Homogeneous  Groupings  (Non-SI  Units 


Run  Plant 


II  0-1  none 


II  0-2  none 


II  0-3  none 


Density 
M.  1/m2 


Mean 
Velocity 
V,  m/sec 


Bed 

Average 

Hydraulic 

Bed 

n 

Radius  m 

V*fV 

0.000  2.355 


0.000 


0.00013 


0.000  14.788  11.940  10.00015  10.016 


1.844 


0.016  2.901 


1.67  0.221  1.690  3.470 


0.00693 


1.300  2.460  0.00496  0.042 


2.33  0.250  3.959  0.963  0.00030  0.042  3.720 


2.33  0.250  3.225  1.570  0.00063  0.035  3.011 


2.33  0.250  3.490  1.934  0.00085  0.034  3.244 


2.33  0.250  3.125 


2.33  0.250 


0.00043  0.045  2.979 


1.699  0.00125  0.040  2.219 


2.33  0.250  2.565  2.013  0.00110  0.033  2.410 


2.33  0.250  2.787  2.270  0.00123  0.032  2.603 


2.33  0.250  2.676  2.522  0.00167  0.033  2.516 


2.33  0.250  2.454  2.827  0.00199  0.031 


3.002  3.102  0.00191  0.030  1 


Purpleleaf 

Euonvmus 


1.190  3.878  1.048  10.00041  10.045 


0.0200 


0.0200 


0.0220 


1  1-1 

Yellow  Twig 
Dogwood 

0.498 

4.170 

1.200 

0.00053 

0.046 

3.944 

0.216 

1  1-2 

Yellow  Twig 
Dogwood 

1.67 

0.498 

4.120 

2.000 

0.00124 

0.042 

3.885 

0.198 

0.0650 

1  1-3 

Yellow  Twig 
Dogwood 

1.67 

0.498 

3.680 

2.460 

0.040 

3.474 

0.185 

0.0590 

1  1-4 

Yellow  Twig 
Dogwood 

1.67 

0.498 

3.090 

1.580 

0.00119 

0.047 

2.959 

0.213 

0.0670 

1  1-5 

Yellow  Twig 
Dogwood 

1.67 

0.498 

3.350 

1.930 

0.00140 

0.043 

3.185 

0.196 

0.0620 

1  1-6 

Yellow  Twig 
Dogwood 

1.67 

0.498 

3.440 

2.260 

0.00163 

0.040 

3.252 

0.183 

0.0580 

1  1-7 

Yellow  Twig 
Dogwood 

1.67 

0.498 

1.760 

2.880 

0.00582 

0.048 

1.710 

0.197 

0.0560 

1  1-8 

Yellow  Twig 
Dogwood 

0.498 

2.350 

3.250 

0.00477 

2.258 

0.181 

0.0540 

1  1-9 

Yellow  Twig 
Dogwood 

1.67 

0.498 

2.910 

3.580 

0.00418 

0.038 

2.766 

0.170 

0.0530 

1  2-1 

Yellow  Twig 
Dogwood 

1.67 

0.221 

4.450 

2.510 

0.00102 

0.031 

4.041 

0.145 

0.0480 

Yellow  Twig 
Dogwood 

1.67 

0.221 

3.770 

3.030 

0.00165 

0.031 

3.463 

0.142 

0.0460 

0.174  0.0500 


0.191  0.0530 


0.195  0.0640 


0.0500 


0.154  0.0490 


0.204  0.0640 


0.176  0.0530 


0.145  0.0440 


0.141  0.0430 


0.146  0.0450 


0.136  0.0410 


0.133 


0.209  | 


1  4-2 

Purpleleaf 

Euonymus 

0.67 

1.190 

3.921 

1.377 

0.00055 

0.040 

3.681 

0.186 

0.0600 

1  4-3 

Purpleleaf 

Euonymus 

0.67 

1.190 

3.673 

2.195 

0.00159 

0.042 

3.489 

0.195 

0.0630 

1  4-4 

Purpleleaf 

Euonymus 

0.67 

1.190 

2.762 

2.172 

0.00225 

0.045 

2.658 

0.202 

0.0620 

1  4-5 

Purpleleaf 

Euonymus 

0.67 

1.190 

2.911 

2.512 

0.00251 

0.042 

2.787 

0.189 

0.0590 

1  4-7 

Purpleleaf 

Euonymus 

0.67 

1.190 

1.610 

2.679 

0.00477 

0.042 

1.565 

0.183 

0.0520 

1  5-1 

Purpleleaf 

Euonymus 

0.67 

0.529 

3.385 

1.348 

0.00053 

0.038 

3.177 

0.172 

1  5-2 

Purpleleaf 

Euonymus 

0.67 

0.529 

3.394 

2.074 

0.00106 

0.035 

3.172 

0.159 

0.0500 

1  5-3 

Purpleleaf 

Euonymus 

0.67 

0.529 

2.320 

3.158 

0.00436 

0.040 

2.231 

0.177 

0.0530 

1  6-1 

3.17 

0.113 

4.143 

1.059 

0.00110 

0.075 

0.357 

0.1190 

1  6-2 

Red  Twig 
Dogwood 

3.17 

0.113 

4.148 

1.573 

0.00213 

0.070 

4.046 

0.336 

0.1110 

mm 

Red  Twig 
Dogwood 

3.17 

0.113 

4.252 

2.005 

0.00266 

0.062 

4.129 

mm 

0.0990 

1  6-4 

3.17 

0.113 

3.085 

1.139 

0.00204 

0.085 

3.036 

0.390 

0.1230 

1  6-5 

Red  Twig 
Dogwood 

3.17 

0.113 

2.485 

1.997 

0.00508 

0.070 

2.442 

0.313 

0.0950 

0.00833  0.070 


0.00540  0.050 


0.00117  0.070 


0.00322  0.070 


0.00145  0.063 


0.00180  0.050 


0.00229  0.043 


0.00276  0.038 


0.00132  0.050 


0.00157  0.042 


0.00019  0.071 


0.00059  0.053 


0.00112  0.040 


0.00201  0.032 


EEMlEEgglil 


0.037 


0.00095  10.035 


0.00103  0.033 


0.00119  0.030 


0.00099  0.083 


0.00163  0.070 


Table  5 

Summary  of  Large  Flume  Results  with  Mixed  Plant  Groupings  (SI  Units) 


Run  Plants 


2-1  20  Service  Berry,  68 

Yellow  Twig  Dogwood, 
68  Euonymus 


20  Service  Berry,  68 
Yellow  Twig  Dogwood, 
68  Euonymus 


2-3  20  Service  Berry,  68 

Yellow  Twig  Dogwood, 
68  Euonymus 


2-4  20  Service  Berry,  68 

Yellow  Twig  Dogwood, 
68  Euonymus 


2-5  20  Sen/ice  Berry,  68 

Yellow  Twig  Dogwood, 
68  Euonymus 


20  Service  Berry,  68 
Yellow  Twig  Dogwood, 
68  Euonymus 


3-1  68  Yellow  Twig 

Dogwood,  68  Euonymu; 


3-4  68  Yellow  Twig  3.66 

Dogwood,  68  Euonymus _ 


2.48 


2.48 


2.48 


lEim 


jg»i  i  mi  1 1  I  'I  I  M 


23  Salt  Cedar,  83  Black  4.20 
Willows,  50  Red  Willows _ 


11-2  23  Salt  Cedar,  83  Black  4.20 
Willows,  50  Red  Willows  _ 


11-3  23  Salt  Cedar,  83  Black  4.20 
Willows,  50  Red  Willows _ 


11-4  23  Salt  Cedar,  83  Black  4.20 
Willows,  50  Red  Willows  _ 


11-5  23  Salt  Cedar,  83  Black  4.20 
Willows,  50  Red  Willows _ 


11-6  23  Salt  Cedar,  83  Black  4.20 
Willows,  50  Red  Willows 


11-7  23  Salt  Cedar,  83  Black  4.20 
Willows,  50  Red  Willows 


83  Black  Willows,  50 
Red  Willows 


Mean 

Velocity  Energy 

V,  m/sec  Slope  S  Average  n 


Table  5  (Concluded 


Run  Plants 


12-2 

83  Black  Willows,  50 

Red  Willows 

3.58 

1.426 

0.551 

0.00113 

0.046 

1.353 

0.220 

0074 

12-3 

83  Black  Willows,  50 

Red  Willows 

3.58 

1.388 

0.763 

0.00210 

0.045 

1 .320 

0.215 

0.072 

83  Black  Willows,  50 

Red  Willows 

3.58 

0.906 

0.910 

0.00333 

0.041 

0.867 

0.186 

0.058 

12-5 

83  Black  Willows,  50 

Red  Willows 

3.58 

0.821 

0.789 

0.00326 

0.045 

0.791 

0.202 

0.062 

12-6 

83  Black  Willows,  50 

Red  Willows 

3.58 

0.776 

0.726 

0.00228 

0.040 

0.743 

0.178 

0.054 

12-7 

83  Black  Willows,  50 

Red  Willows 

3.58 

0.680 

0.688 

0.00175 

0.035 

0.647 

0.151 

0.045 

Table  6 

Summary  of  Large  Flume  Results  with  Mixed  Plant  Groupings  (Non-SI  Units) 


2-1  20  Service  Berry,  68 

Yellow  Twig  Dogwood, 
68  Euonymus 


20  Service  Berry,  68 
Yellow  Twig  Dogwood, 
68  Euonymus 


2-3  20  Service  Berry,  68 

Yellow  Twig  Dogwood, 

68  Euonymus  _ 


2-4  20  Service  Berry,  68 

Yellow  Twig  Dogwood, 
68  Euonymus 


2-5  20  Service  Berry,  68 

Yellow  Twig  Dogwood, 
68  Euonymus 


20  Service  Berry,  68 
Yellow  Twig  Dogwood, 
68  Euonymus 


3-1  68  Yellow  Twig 

Dogwood,  68  Euonymus 


68  Yellow  Twig 
Dogwood,  68  Euonymus 


11-1  23  Salt  Cedar,  83  Black  0.39 
Willows,  50  Red  Willows _ 


11-2  23  Salt  Cedar,  83  Black  0.39 
Willows,  50  Red  Willows _ 


11-3  23  Salt  Cedar,  83  Black  0.39 
Willows,  50  Red  Willows _ 


11-4  23  Salt  Cedar,  83  Black  0.39 
Willows,  50  Red  Willows  _ 


11-5  23  Salt  Cedar,  83  Black  0.39 
Willows,  50  Red  Willows _ 


11-6  23  Salt  Cedar,  83  Black  0.39 
Willows,  50  Red  Willows _ 


11-7  23  Salt  Cedar,  83  Black  0.39 
Willows,  50  Red  Willows 


12-1  83  Black  Willows,  50 
Red  Willows 


Table  7 

Small  Flume  Results,  Drag 


Plant 


Staahom  Sumac 


Staahom  Sumac 


Staahom  Sumac 


Staahom  Sumac 


Staahom  Sumac 


Staahom  Sumac 


Staahom  Sumac 


Staahom  Sumac 


Staahom  Sumac 


Staahom  Sumac 


Artie  Blue  Willow 


Artie  Blue  Willow 


Artie  Blue  Willow 


Artie  Blue  Willow 


Artie  Blue  Willow 


Artie  Blue  Willow 


Artie  Blue  Willow 


Artie  Blue  Willow 


Artie  Blue  Willow 


Artie  Blue  Willow 


Norway  Maple 


Norway  Maple 


Norway  Maple 


Norway  Maple 


Norway  Maple 


Norway  Maple 


Norway  Maple 


Norway  Maple 


Norway  Maple 


Norway  Maple 


Western  Sand  Che 


Western  Sand  Cher 


Western  Sand  Cher 


Western  Sand  Cher 


Western  Sand  Che 


Western  Sand  Che 


Western  Sand  Cheri 


Western  Sand  Cheri 


Western  Sand  Cher 


Western  Sand  Cher 


Common  Privet 


Common  Privet 


Common  Privet 


Common  Privet 


Common  Privet 


Common  Privet 


Common  Privet 


Common  Privet 


Common  Privet 


Common  Privet 


Measurements _ 


With  Leaves 


roach  Velocl 


ft/sec  m/sec 


1.63  0.50 


2.15  0.66 


am 


Blue  Elberber 


Blue  Elberbe 


Blue  Elberberi 


Blue  Elberber 


Blue  Elberber 


Table  7  (Concluded 


2S9 

Jul 

m 

m 


Plant 


Blue  Elberbe 


Blue  Elberbe 


Blue  Elberber 


Blue  Elberber 


French  Pink  Pussywillow 


French  Pink  Pussywillow 


French  Pink  Pussywillow 


French  Pink  Pussywillow 


French  Pink  Pussywillow 


French  Pink  Pussywillow 


French  Pink  Pussywillow 


With  Leaves 


roach  Veloc 


m/sec  Ibf 


1.03  2.406 


1.19  12.783 


Without  Leaves 


roach  Velocl 


Sycamore 


Sycamore 


Sycamore 


Sycamore 


Sycamore 


Sycamore 


Sycamore 


Yellow  Twiq  Doqwood 


Yellow  Twiq  Doawood 


Yellow  Twiq  Doqwood 


Yellow  Twiq  Doqwood 


IKmMGEnESJE 


Yellow  Twiq 


Yellow  Twiq 


Yellow  Twiq 


Yellow  Twiq 


Yellow  Twiq 


Yellow  Twiq 


Yellow  Twiq 


Yellow  Twiq 


Yellow  Twiq 


Yellow  Twiq 


Doqwood 


Doqwood 


Doqwood 


Doqwood 


Doqwood 


Doqwood 


Doqwood 


Doqwood 


Doqwood 


Doqwood 


373 

373 

373 

373 

373 


Yellow  Twiq  Doqwood 


Yellow  Twiq  Doqwood 


Yellow  Twiq  Doqwood 


Yellow  Twiq  Doqwood 


Euonvmus 


Euonvmus 


Euonvmus 


Euonvmus 


Yellow  Twiq  Doqwood 


Yellow  Twiq  Doqwood 


Yellow  Twiq  Doqwood 


Yellow  Twiq  Doqwood 


Yellow  Twia  Doawood 


IgEmEEEjggSETI 


Force 


N 


3.189 


4.555 


6.379 


8.856 


0.854 


1.281 


1.668 


2.011 


2.353 


3.723 


4.493 


0.258 


0.427 


0.601 


0.814 


1.028 


1.090 


1.219 


2.011 


2.353 


2.460 


0.480 


0.916 


1.308 


1.833 


2.006 


2.006 


2.771 


2.789 


2.922 


2.616 


Table  8 

Large  Flume  Results,  Drag 


Measurements 


[33 

m: 


5 


6 


7-1:1 


7-1 :2 


3 


7-1:4 


8-1:1 


8-1:2 


8-1:3 


8-1:4 


8-1:5 


8-1 :6 


8-1:7 


11-1:1 


2 


11-1:3 


11-1:4 


11-1:5 


11-1:6 


11-1:7 


11-1:8 


11-1:9 


11-1:10 


Plant 


Mulefat 


Mulefat 


Mulefat 


Mulefat 


Mulefat 


Mulefat 


Alder 


Alder 


Alder 


Alder 


Valiev  Elderbe 


Valiev  Elderberi 


Valiev  Elderberi 


Valiev  Elderbe 


Black  Willow 


Black  Willow 


Black  Willow 


Black  Willow 


Black  Willow 


Black  Willow 


Black  Willow 


Black  Willow 


Black  Willow 


Black  Willow 


i  Leaves 


roach  Velocl 


Table  9 

Modulus  of  Plant  Stiffness  for  Evaluated  Plants 


Plant  name 


Alder 


Arctic  Blue  Willow 


Black  Willow 


Blue  Elderbe 


Common  Privet 


Yellow  Twig  Dogwood 


Red-osier  Dogwood 


Berried  Elderber 


Puroleleaf  Euonvmus 


Mountain  Black  Willow 


Mulefat 


Norwav  Maple 


French  Pink  Pussywillow 


Red  Willow 


Salt  Cedar 


Sen/ice  Be 


Staghorn  Sumac 


Sycamore 


Valiev  Elderberi 


Western  Sand  Cher 


Sand  Bar  Willow 


Pacific  Willow 


Lemon's  Willow 


Wild  Rose  Bush 


iemieu 


Salix  purpurea  nana 


Salix  nigra 


Sambucus  Canadensis 


Ligustrum  vuiaare 


Comus  stolonifera  flaviramea 


Comus  Sericea 


Sambucus  racemosa 


Euonvmus  fortune i  colorata 


Saiix  monticoia 


Baccharis  alutinosa 


Acer  platenoides 


Salix  caprea  pendula 


Salix  s, 


Tamarix  s 


Amelanchier 


Rhus  tvphina 


Platenus  acer  if olia 


Sambucus  mexicana 


Prunis  besse 


I  Salix  exigua 


Salix  lasiandra 


Salix  lemonii 


Rosas 


1.804e+06 


4.091  e+05 


2.930e+05 


7.7040e+05 


2.550e+06 


4.342e+06 


8.168e+04 


2.278e+06 


7.430e+05 


8.992e+05 


4.569e+06 


3.345e+05 


8.81 0e+05 


3.930e+06 


4.003e+06 


1 .095e+06 


3.244e+06 


7.672e+06 


3.567e+06 


4.990e+06 


5.300e+06 


4.090e+06 


6.070e+06 


3.768e+04 


8.544e+03 


6.119e+03 


5.708e+03 


1 .609e+04 


5.326e+04 


9.069e+04 


1.706e+03 


4.758e+04 


1.552e+04 


IKE 2E35D 


9.542e+04 


6.986e+03 


1.840e+04 


8.207e+04 


8.360e+04 


2.288e+04 


6.774e+04 


1.602e+05 


7.449e+04 


1.120e+05 


8.530e+04 


Appendix  A 
Example  Problem 


The  equations  developed  in  this  report  allow  hydraulic  roughness  values  to 
be  determined  for  homogenous  and  nonhomogeneous  flood  plains.  When  a 
number  of  species  are  present  in  the  floodplain  or  area  of  interest,  the  values  for 
the  various  plants  are  either  combined,  as  shown  in  Equation  29,  or  the  flood 
plain  is  broken  into  homogenous  areas  that  are  then  either  solved  simultaneously 
for  the  flow  or  aggregated  to  provide  a  representative  roughness  (i.e.. 

Manning’s  n)  value. 

The  following  is  an  example  for  determining  the  vegetative  resistance  and 
the  equivalent  resistance  of  the  left  bank  of  a  floodplain.  The  left  bank  is 
divided  into  three  subareas  with  the  far  left  area  vegetated  with  shrubs,  the 
middle  area  vegetated  with  three  different  plant  types,  and  the  right  subarea 
vegetated  with  willows  as  shown  in  Figure  Al.  For  a  given  flow  depth  in  the 
main  channel.  Table  Al  shows  the  flow  depth,  area,  wetted  perimeter,  and 
hydraulic  radius  for  each  subarea.  The  slope  of  each  area  is  0.0002,  the  fluid 
density  is  1,000  kg/m3,  and  the  kinematic  viscosity  is  1.3E+06  m2/s.  The 
properties  of  the  channel  are  summarized  in  Table  Al  for  ease  of  reference. 


Figure  Al .  Left  bank  floodplain  for  example  problem  showing  calculation 
subareas  based  on  vegetation  types 


Appendix  A  Example  Problem 


Al 


Table  Al  | 

Channel  Properties  for  Example  Problem  | 

|[  Channel  Properties 

Subarea  1 

Subarea  2 

Subarea  3 

Overbank  Total 

S 

0.0002 

0.0002 

0.0002 

.0002 

Y0  m 

0.61 

1.524 

2.134 

Flow  Area 

27.87 

46.45 

32.52 

106.84 

P  m 

31.39 

31.09 

15.85 

78.33 

Rh  m 

0.89 

1.49 

2.05 

1.364 

Table  A2  lists  the  plant  characteristics  and  properties  for  the  shrubs  in 
subareas  1  through  3.  The  stiffness  modulus  was  measured  in  the  field  by 
applying  Equation  33  for  subarea  1,  by  using  Equation  35  for  weighted  plant 
characteristic  values  for  subarea  2,  and  by  using  the  value  for  Pacific  willows 
from  Table  9  for  subarea  3  multiplied  by  the  (H/Ds)1 5  value  for  the  willows  in 
subarea  3.  The  use  of  the  value  from  Table  9  and  the  H/D  parameter  for  the 
willows  in  subarea  3  provides  an  estimate  of  the  actual  Es  for  the  willows  in  the 
subarea  that  are  different  in  size  from  those  evaluated  in  the  compilation  of  the 
data  presented  in  Table  9. 


Table  A2  also  lists  the  plant  characteristics  and  shows  the  calculations  to 
determine  the  weighted  average  characteristics  to  be  used  for  the  resistance 
calculations.  This  is  done  by  multiplying  the  various  plant  characteristics  by  that 
plant’s  relative  density  on  the  floodplain  in  decimal  form,  HiMi/Mtotai,  where  the 
subscript  1  refers  to  plant  type  1.  These  values  are  then  summed  for  all  three 
plant  types  to  arrive  at  the  weighted  value  for  H.  Thus,  the  equation  for  the 
weighted  value  becomes  as  follows: 


H 


ave 


Hi 


M  total 


+ 


H  2 


Mi 

M  total 


+ 


H  3 


Mi 

M  total 


This  averaging  technique  was  verified  using  flume  data  from  the  Utah  Water 
Research  Laboratory  Study  and  provides  the  correct  average  values  for  use  in  the 
equations  presented  herein. 

Equation  23  was  used  to  calculate  V*/V  for  the  partially  submerged 
vegetation  in  subarea  1  while  Equation  21  was  used  for  the  fully  submerged 
subareas  2  and  3.  Manning’s  n  values  were  calculated  in  the  same  way  using 
Equation  24  for  subarea  1  and  Equation  22  for  subareas  2  and  3. 


The  calculated  Manning’s  n  values  for  subareas  1  through  3  are  0.075, 
0.088,  and  0.127,  respectively,  while  the  V./V  values  are  0.240,  0.257,  and 
0.353,  respectively.  The  values  for  the  various  subareas  can  then  be  used 
individually  in  a  hydraulic  model  or  composited  to  obtain  a  value  for  the  entire 
channel.  Methods  described  in  the  SAM  User’s  Manual  (Thomas  et  al.  in 
preparation)1  are  recommended  for  compositing  the  subareas  if  that  is  necessary. 


1  References  in  this  appendix  are  cited  in  the  References  section  at  the  end  of  the  main 
text. 
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Table  A2 

Plant  Characteristics  and  Roughness  Calculations  for  Subareas  1-3 

Plant  Characteristics 

Sub  area  1 

Subarea  2 

Subarea  3 

(Pacific 

Willow) 

Plant 

Type  1 

If 

mwi 

Weighted 

Values 

H  m 

0.91 

0.2 

0.71 

0.99 

0.59 

1.83 

hP  m 

0.76 

0.2 

0.51 

0.81 

0.49 

1.52 

We  m 

0.76 

0.254 

0.18 

0.86 

0.48 

1.22 

No.  of  Stems 

4 

2 

6 

1 

2.20 

4 

Ds  cm 

1.27 

0.64 

0.64 

1.2 

0.86 

1.3 

As  m  n  Ds2  /  4  x  (no.  of  stems) 

5.06E-04 

6.43E-05 

1.93E-04 

1.13E-04 

1.03E-04 

5.31  E-04 

0.52 

1.83 

0.61 

1.62 

4.06 

0.11 

M  /  M  total  (%) 

45.0 

15.0 

39.9 

100.0 

H/Ds 

72 

31 

111 

82 

64 

141 

F45.N 

8.08 

Es,  N/m^,  Eq  33 

1.75E+09 

Es.  N/rrf,  Eq  35 

1 .33E+09 

3.15E+08 

3.59E+09 

1.49E+09 

1.28E+09 

6.56E+09 

I  5  f/TiUh 

8.85E+09 

Vegetative  Resistance 

Yo/H 

0.67 

2.72 

1.17 

Submerqence  P  <  0.8H  <  Full 

Partial 

Full 

Full 

Aj  or  A*  (m^) 

0.35 

0.24 

1.85 

V*  =  (qRhS) 

4.18E-02 

5.41  E-02 

6.34E-02 

1449447 

207582 

630112 

MAs* 

0.182 

0.964 

0.204 

V*  Rh  /  v 

28608 

61971 

100009 

WV  (Eq  23) 

.130 

WV  (Eq  21) 

0.380 

0.350 

Manning’s  n  (Eq  22  or  24) 

IISEZni 

0.126 
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